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Abstract 

This paper is concerned with the fading multiple-input multiple-output multiple-access channel (MIMO- 
MAC) with multiple receive antennas at the base station (BS) and multiple transmit antennas at each mobile 
terminal (MT). Two multiple-access techniques are considered for scheduling transmissions from each MT to the 
BS at the same frequency, which are space -division multiple-access (SDMA) and time-division multiple-access 
(TDMA). For SDMA, all MTs transmit simultaneously to the BS and their individual signals are resolved at the 
BS via multiple receive antennas while for TDMA, each MT transmits independently to the BS during mutually 
orthogonal time slots. It is assumed that the channel-state information (CSI) of the fading channel from each 
MT to the BS is unknown at each MT transmitter, but is perfectly known at the BS receiver. Thereby, the BS can 
acquire the long-term channel-distribution information (CDI) for each MT. This paper extends the well-known 
transmit-covariance feedback scheme for the point-to-point fading MIMO channel to the fading MIMO-MAC, 
whereby the BS jointly optimizes the transmit signal covariance matrices for all MTs based on their CDI, and 
then sends each transmit covariance matrix back to the corresponding MT via a feedback channel. The main 
goal of this paper is to characterize the so-called multiuser power region under the multiuser transmit-covariance 
feedback scheme for both SDMA and TDMA. The power region is defined as the constitution of all user transmit 
power-tuples that can achieve reliable transmissions for a given set of user target rates. Simulation results show 
that SDMA can achieve substantial power savings over TDMA for the fading MIMO-MAC, even when the 
number of antennas at the BS is equal to that at each MT. Moreover, this paper demonstrates the usefulness of 
the multiuser power region for maintaining proportionally-fair power consumption among the MTs. 
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I. Introduction 

Transmission through multiple transmit and multiple receive antennas, or the so-called multiple- 
input multiple-output (MIMO) technology, is known as an efficient means for providing enormous 
information rates in rich-scattering mobile environments [l]-[3]. Characterization of the fading MIMO 
channel capacity limits, under various assumptions on the transmitter-side and receiver-side channel-state 
information (CSI) and channel-distribution information (CDI), has motivated a great deal of valuable 
scholarly work (e.g., [4] and references therein). In particular, the case where the CSI is perfectly known 
at the receiver but unknown at the transmitter has drawn much interest due to its validity in many 
practical situations. This is because the presumption of perfect receiver-side CSI is usually reasonable 
for wireless channels where the receiver can locally estimate the fading channel, while the complete 
CSI feedback from the receiver to the transmitter is difficult or even impossible. Consequently, many 
schemes that exploit various forms of partial channel feedback have been proposed in literature. Among 
others, the transmit-covariance feedback scheme is known to be capable of achieving data rates close to 
the fading MIMO channel ergodic capacity when the channel CDI exhibits some long-term consistent 
statistical properties, e.g., constant channel mean and/or constant channel covariance matrix [5]-[8]. In 
the transmit-covariance feedback scheme, the receiver determines the transmit signal covariance matrix 
based on the CDI, and then sends it back to the transmitter through a feedback channel. In [5]-[8], 
the feedback transmit covariance matrix for optimizing the channel ergodic capacity, and the conditions 
under which the beamforming - the transmit covariance matrix has rank one - is optimal, have been 
established for the single-user fading channel. In this scheme, the transmit covariance matrix is fixed 
as long as the CDI is not changed. Therefore, this scheme requires much less feedback complexity and 
is also more robust to the delay of the feedback channel compared to other partial channel feedback 
schemes based on the instantaneous MIMO channel realizations (e.g., [9], [10] and references therein). 

This paper considers the fading MAC with additive white Gaussian noise (AWGN) at the receiver, 
and assumes that the CSI from each mobile terminal (MT) to the base station (BS) is unknown at each 
MT transmitter, but is perfectly known at the BS receiver. Thus, the BS can acquire the channel CDI for 
each MT. This paper extends the transmit-covariance feedback scheme for the single-user fading MIMO 
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channel to the fading MIMO multiple-access channel (MIMO-MAC) where multiple antennas are used 
by the BS and possibly by each MT. Two multiple-access techniques are considered for scheduling 
transmissions from each MT to the BS at the same frequency: space-division multiple-access (SDMA) 
and time-division multiple-access (TDMA). For SDMA, all MTs transmit simultaneously to the BS and 
their individual signals are decoded jointly at the BS while for TDMA, each MT transmits independently 
to the BS during mutually orthogonal time slots and thus only single-user decoding is needed. The 
multiuser transmit-covariance feedback scheme is then described as follows. For SDMA, the BS first 
jointly optimizes the transmit signal covariance matrices for all MTs, based on the multiuser CDI as 
well as the rate requirement and the power budget of each MT, and then sends them back to each MT 
for transmission. This scheme has also been considered in [11], [12] for characterizing the capacity 
region and establishing the conditions for the optimality of beamforming for the fading MIMO-MAC, 
respectively. In contrast, for TDMA, the BS jointly optimizes the duration of transmission time slot for 
each MT along with their transmit covariance matrices. These optimized values are then sent back to 
each corresponding MT via the feedback channel. 

This paper studies the information-theoretic limits of the fading MIMO-MAC under the multiuser 
transmit-covariance feedback scheme when either SDMA or TDMA is employed. Two commonly 
adopted means to measure the information-theoretic limits of multiuser channels are the capacity region 
and the power region. The capacity region is defined as the constitution of all achievable rate-tuples 
for the users given their individual power constraints, while the power region consists of all possible 
power-tuples for the users under which a given rate-tuple is achievable. This paper is mainly concerned 
with the characterization of the multiuser power region. Our motivations are justified as follows: 

First, characterization of the power region for the fading MIMO-MAC is a challenging problem. 
Considering initially the case of SDMA, the capacity region of a deterministic (no fading) Gaussian 
MAC with a single transmit and a single receive antenna (SISO-MAC) has the well-known polymatroid 
structure [13], which also holds for the fading MIMO-MAC. On the other hand, the power region of 
a deterministic SISO-MAC is known to have a contra-polymatroid structure [14]. The polymatroid and 
the contra-polymatroid structures can be utilized to reduce significantly the computational complexity of 
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finding the boundary points of the capacity region and the power region, respectively [13], [14]. However, 
the contra-polymatroid structure is non-existent for the power region when the channel exhibits fading 
[15] and/or the BS uses multiple antennas [16]Q As a result, characterization of the power region for the 
fading MIMO-MAC under SDMA is yet fully understood in literature. On the other hand, for TDMA, 
given the duration of each MT transmission time slot (e.g., equal time-slot durations for all MTs in 
the conventional TDMA), the BS only needs to optimize the transmit covariance matrices for the MTs 
independently such that their individual transmission powers are minimized for supporting their own 
target rates. However, with time-slot duration adjustment for each MT, the BS now needs to consider 
the more challenging problem of jointly optimizing the time-slot durations and the transmit covariance 
matrices for all MTs. This joint optimization has been less studied in literature. 

Secondly, characterization of the power region can potentially provide very useful insights on the 
resource allocation problem for wireless networks, e.g., the wireless cellular network. In many practical 
situations, each BS in the cellular network controls the transmit power of each MT in its cell such 
that each individual MT rate demand - transmission quality-of-service (QoS) - is satisfied (e.g., [17], 
[18]). Power control can be beneficial in many aspects, e.g., to maintain fair power consumption among 
MTs, to tailor for each MT's peak-power constraint, and to mitigate the co-channel interference between 
multiple cells operating at the same frequency. By exploiting the multiuser power region, the minimum 
power consumption in the network can be achieved under practical transmission constraints. 

The main contributions of this paper are summarized as follows: 

• For both SDMA and TDMA, the paper presents efficient algorithms for characterizing each bound- 
ary point of the power region for the fading MIMO-MAC. The developed algorithms are based on a 
Lagrange primal-dual approach that exploits a novel dual relationship between the power region and 
the corresponding capacity region for the fading MIMO-MAC. For SDMA, the proposed algorithm 
determines jointly the optimal transmit covariance matrices for all MTs as well as their optimal 
decoding order at the receiver. For TDMA, all MT transmit covariance matrices along with their 
assigned variable time-slot durations are jointly optimized. 

'More discussions on this aspect are postponed to Section HV- Al 
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• In addition to the conventional way to characterize the boundary of the power region by solving a 
sequence of user weighted sum-power minimization problems subject to fixed user rate constraints, 
this paper proposes an alternative means for such characterization by considering the user sum- 
power minimization problem under different user power-profile constraints, where the user power- 
profile regulates the power consumption of users under some given proportional fairness. 
The remainder of this paper is organized as follows. Section [II] introduces the system model for the 
fading MEVIO-MAC and describes the proposed multiuser transmit-covariance feedback scheme under 
SDMA and TDMA. Section [III] provides the definition of the power region for the fading MIMO- 
MAC, together with two problem formulations for characterization of the power region, based on the 
user weighted sum-power minimization and the user power-profile vector, respectively. Section [IV] and 
Section [V] study the power region in the case of SDMA and TDMA, respectively, and present efficient 
algorithms for characterizing the power region in each case. Section [VI] provides numerical results to 
verify the usefulness of the proposed algorithms under realistic channel parameters. Finally, Section |VII] 
concludes the paper. 

Notations: Scalar signals are denoted by lower-case letters, e.g., x, y. Bold-face lower-case letters are 
used for vector signals, e.g., x, y, and bold-face upper-case letters for matrices, e.g., S, M. \S\ denotes 
the determinant, S~ x the inverse and Tr(S') the trace of a square matrix S. For any general matrix M, 
M T and denote its transpose and conjugate transpose, respectively, i" denotes the identity matrix. 
E[-] denotes statistical expectation. C xxy denotes the space of x x y matrices with complex entries. R M 
denotes the M -dimensional real Euclidean space and is the nonnegative orthant. The distribution 
of a circular symmetric complex Gaussian (CSCG) vector with the mean vector x, and the covariance 
matrix £ is denoted by CJ\f(x, £), and ~ means "distributed as." The sign >z denotes the generalized 
inequality [19] and for a square matrix S, S >z means that S is positive semi-definite. mm(x,y) 
denotes the minimum between two real numbers x and y. 

II. System Model 

This paper considers a narrow-band flat fading MIMO-MAC with r receive antennas at the BS and 
K MTs equipped with ti, . . . ,t K transmit antennas, respectively. All MTs transmit synchronously to 
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the BS by sharing a common frequency band. It is assumed that the space of fading states is continuous 
and infinite, and the fading process is stationary and ergodic. Under the assumption that the transmitted 
symbol period is equal to the inverse of the common transmission bandwidth for all MTs, at each fading 
state v, the equivalent discrete-time MAC is given by 

K 

y = s ^H k (v)x k + z, (1) 

fe=i 

where y E C rxl denotes the received signal vector, x k E C ifcXl and H k (u) E C rxtk denote, respectively, 
the transmitted signal vector and the channel matrix of MT k, k = 1,...,K; z E C rxl denotes the 
vector of additive noise at the receiver, and it is assumed that z ~ CJ\f(0, 1). 

This channel model also assumes that the CSI is perfectly known at the BS but is unknown at each 
MT. With the CSI available, the BS can acquire the long-term CSI statistics (or equivalently, the CDI) of 
each MT. Based on the multiuser CDI, the BS determines the transmit signal covariance matrices for all 
MTs jointly according to their individual rate requirement and power budget, and then sends them back 
to each MT for transmission. This paper refers to this scheme as multiuser transmit-covariance feedback. 
Let the transmit covariance matrix of MT k be S k — ^\x k x k ], where the expectation is taken over the 
code-book and S k >z 0. S k is assumed to be fixed for all fading states v. Fig. \T\ illustrates the system 
model considered in this paper. Since this paper is concerned with the information-theoretic limits of 
a Gaussian MAC, the optimal Gaussian code-book is assumed for each MT, i.e., x k ~ CJ\f(0, S k ),Wk. 
The transmit covariance matrix of MT k can be expressed by its eigenvalue decomposition as 

S k = V k ^ k V[. (2) 

V k E C tkXdk is known as the precoding matrix where V\V k = I, and d k < min(t fc ,r). d k is usually 
referred to as the spatial multiplexing gain as it measures the number of degrees of transmission freedom 
(or equivalently, the number of data streams) in the spatial domain. If d k is equal to one, the associated 
transmission scheme is usually referred to as beamforming. Yi k is a d k x d k diagonal matrix with positive 
diagonal elements that provide the power-loading to the corresponding transmitted data streams. The 
transmitter of each MT can be implemented as the cascade of the following operations: encoding the 
information bits by the optimal Gaussian code-book, interleaving coded symbols randomly into each data 



C smk ({S k }) = \ r E R£ : J2 r k < ^ 



~ log 



VJC{1,...,K} . (3) 



stream, power-loading and then jointly precoding all data streams. Next, two multiple-access techniques 
considered in this paper are illustrated, namely, SDMA and TDMA. 

For SDMA, all MTs transmit simultaneously to the BS. In this paper, it is assumed that the BS receiver 
uses the optimal (capacity-achieving) multiuser detection. For a fixed set of {S k }, k = 1, . . . , K, all 
the rate-tuples in the set, Csdma defined below, are achievable (e.g., [11], [12]): 

Y,H k {v)S k H\{v)+I 

The code-book of MT k should satisfy Tr(S k ) = Pk, where p = (pi, . . . ,px) € denotes the vector 
of average transmit powers for the MTs. 

On the other hand, for TDMA, the BS divides the total transmission time into multiple transmission 
blocks of equal duration T. Each transmission block is then further divided into K non-overlapping 
time slots assigned to the MTs. These time slots are assumed fixed over all blocks. During the time slot 
of MT k, only this MT communicates with the BS and other MTs are silent, i.e., x k ' = 0, VA;' ^ k in 
(OQ). Let r k T denote the time-slot duration assigned to MT k, where < r k < l,Wk and J2 k =i T k — !• 
The BS determines jointly the slot duration for each MT and their transmit covariance matrices, and 
then sends them back to each MT. For a fixed set of {r k } and {S k }, k = 1, . . . , K, each MT transmits 
over a single-user fading MIMO channel studied in e.g., [1], [2], and thus the following rate-tuples in 
the set, C T unk({Tk}, {Sk}), defined below are achievable: 



C Tmk ({r k }, {S k }) = jr G R£ : r k < r k E u 



- log 

2 & 



H k (u)^Hl(u)+I 



\fke{l,...,K}\. (4) 



And, again, p = (p 1 , . . . ,p K ) E IR^ denotes the average transmit powers for the MTs. Notice that for 
TDMA, the actual transmit power for MT k during its assigned time-slot duration r k T is Tr (S k ) /r k , 
but the average transmit power p k over each block duration T is Tr(S k ), the same as SDMA. 

Remark 2.1: In this work, for both SDMA and TDMA, we consider each user's achievable rate in 
the "ergodic" sense, i.e., averaged over all ergodic fading states. In the case of fast-fading channel, the 
resultant ergodic capacity can be achievable by assigning each MT a constant-rate code-book for which 
the codeword length is sufficiently long so as to exploit the ergodicity of the channel. In contrast, in the 
case of slow-fading channel, each MT's codeword might not be able to span over all possible fading 
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states because of practical transmission delay constraint. However, if each MT uses multiple code-books 
with variable rates, the BS, based on the instantaneous channel, can determine the transmission rate 
of each MT and then sends back the corresponding code-book index to each MT for transmission. As 
in the fast-fading case, the same ergodic capacity (sometimes known as the expected capacity) can be 
achieved for each MT via time-averaging its transmission rates over different fading states. 

III. Power Region for Fading MIMO-MAC: Definitions and Characterizations 

The multiuser power region for the fading MAC is defined as the constitution of all user power-tuples 
under which a given set of rates is achievable for all the MTs. Let R = . . . , Rk) G denote 
the vector of rate requirements for the MTs. The power region is then defined as 



Psdma (R) = {pER«:3 {S k }, such that R G C smA ({S k }) ,p k = Tr(S k ), k = l,...,K}, (5) 
for SDMA, and 

^tdma(^) = {p G : 3 {r k }, {S fe },such that R G C TDMA ({r fc }, {S k }) ,p k = Tr(S k ),k = 1, . . . ,K) , (6) 



for TDMA. It is not hard to show that the power regions for both SDMA and TDMA are convex sets. 
The power region is illustrated in Fig. |2] for a two-user fading MAC under either SDMA or TDMA. The 
solid line in Fig. |2] represents the boundary points of the power region, which correspond to all pareto 
optimal power-tuples each minimizing a weighted-sum of the powers among the MTs. Each boundary 
point, e.g., point A as indicated in Fig. [2J might be characterized by two alternative means described 
as follows. 

First, because of the convexity of the power region, each boundary point can be expressed as the 
solution to a weighted sum-power minimization (W-SPmin) problem stated below, for some nonnegative 
user weights, X k , k = 1, . . . , K. For SDMA, the W-SPmin problem can be expressed as 



Problem 1: 



K 



{r k },{S k } 

Subject to 



Minimize 




(7) 



fc=i 



r k > R k \/k 



(8) 



r e C S mk({S k }) 



(9) 



s k y o Vfc. 



(10) 
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Notice that {r k } are auxiliary variables. For TDMA, this W-SPmin problem is given by 



Problem 2: 



K 



{r k },{S k },{r k } 
Subject to 



Minimize 



^A fc Tr(Sy 



(11) 



k=l 



rk > Rk VA; 



(12) 



r e C Tmk ({r k }, {S k }) 



(13) 



r k > VA; 



(14) 



K 



5> = i 



(15) 



k=l 



s k yo vfc. 



(16) 



By definition of the power region for a given target rate-tuple R, each boundary point of the power 
region, p*, under either SDMA or TDMA, can be expressed as a power-tuple supporting the target rates 
in R that has the minimum weighted-sum, Ylk=\ ^kP* k > f° r some weight vector A, among all the power- 
tuples that can support R. Alternatively, the connection between the power-tuple p* on the boundary 
of the power region and the target rate -R can be established by employing the capacity region concept. 
Fig- 12] gives an illustration for this important observation. Consider the power region of a two-user 
MAC under the rate constraint (R 1 ,R 2 ), as shown in Fig. [3] (a). Given Ai and A 2 , the solution to the 
W-SPmin problem is denoted by (pl,p 2 ), which is represented by point A in Fig. |3] (a), and satisfies 
A1P1+A2P2 = P*, where p* is the minimum value of the W-SPmin problem. Thus, the required minimum 
power-pair is (Pi,p 2 ) f° r -^2)- On the other hand, it is not hard to verify that the rate-pair (Ri, R 2 ) 
is on the boundary of the capacity region for the same MAC under the weighted sum-power constraint 
given by A1P1 + \ 2 p 2 < P* ■ This is shown by point B in Fig. [3](b). Moreover, because of the convexity 
of the capacity region, R 2 ) must be the solution to a weighted sum-rate maximization (W-SRmax) 
problem for a given nonnegative user weight vector p, as shown in Fig. [3](b). The above observation 
has an important consequence, i.e., each power-region boundary point as the solution to the W-SPmin 
problem can be equivalently characterized as a boundary point of the corresponding capacity region 
under the same weighted sum-power constraint. As will be shown later in this paper, this result also 
motivates the proposed algorithms for the W-SPmin problem under both SDMA and TDMA. 

Alternatively, each boundary point of the power region can also be considered geometrically as the 
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intersection of a line passing through the origin (the power-tuple with all zeros) and the boundary of 
the power region (see point A in Fig. [2]). Let each line be characterized as p k = a k P, for k = 1, . . . , K 
and P > 0. The vector ex. = (a 1 , . . . , ax) E IR+ is referred to as the user power-profile vector in this 
paper, and it is assumed that J2k=i a fe = 1- The point where the line specified by a intersects the power 
region boundary can be then obtained as the solution to the following optimization problem, referred 
to as the sum-power minimization under the power-profile constraint (SPmin-PPC). For SDMA, this 
problem can be expressed as 
Problem 3: 

Minimize P (17) 

P,{r k },{S k } 

Subject to r k > R k Wk (18) 

r e C smA ({S k }) (19) 

S k y Wk (20) 

Tr(S fc ) <a k P Vk (21) 

P > 0. (22) 
And similarly, the SPmin-PPC problem can be defined for TDMA. The characterization of the power 
region via some prescribed power-profile vector is useful for the BS to regulate the power consumption 
of MTs in a desired proportionally-fair manner. 

It is not hard to show that both problems, W-SPmin and SPmin-PPC, are convex for either SDMA 
or TDMA, and hence, they can be solved by applying convex optimization techniques. The following 
two sections present the solutions to these problems for SDMA and TDMA, respectively. 

IV. Characterization of Power Region for SDMA 
A. Solutions to Weighted Sum-Power Minimization 

This subsection considers Problem[TJ the W-SPmin problem in the case of SDMA. For the special case 
of a deterministic (no fading) SISO-MAC, the W-SPmin can be simplified using the contra-polymatroid 
structure as proposed in [13], [14]. However, as shown later in this subsection, the approach in [13], 
[14] can not be applied here to handle the more general case of the fading MIMO-MAC. Thus, an 
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alterative approach is proposed. 

Deterministic SISO-MAC: Consider a deterministic SISO-MAC where the channel gain for each 
MT in ([T]) is a positive constant, i.e., \H k (v)\ = \fh~ k ^v. In this case, the solutions to the W-SPmin 
problem are obtained as follows. 

Theorem 1: For a deterministic SISO-MAC consisting of K users with channel gains hi, ... , hx, 
and rate requirements Ri, . . . , Rk, the solutions to the W-SPmin problem under SDMA are given by 
[14, Lemma 3.2]: 

f exp ( 2 ^a))-i if Jfc = l 

Pn ^ I exp(2 Et i A, (i) ) -ex P (2 ) ^ 

I ^(fe) ' ' ' ' ' ' 

where the permutation % indicates the optimal decoding order (user ir(l) is decoded last and user %(K) 
is decoded first) at the receiver according to 

Vi>...>^EW. (24) 
Fig. ID shows the connection between each boundary point on the power region and the corresponding 
capacity region earlier described in Section [Till for a two-user deterministic SISO-MAC. Fig. [4] (a) shows 
the power region given rate constraint (R^R?). Consider a vertex (pi, pi) of this power region, and 
arbitrary positive weights Ai, A 2 such that ^ > According to Theorem [H (??*, * s me optimal 
solution to the W-SPmin problem, and is achievable by decoding order 2 — > 1 (user 2's message is 
decoded before user l's). For this given choice of Ai and A 2 , the boundary curve of the capacity region 
under the weighted-sum power constraint AiPi + A 2 £>2 < P*, shown in Fig. |4] (b), can be represented 
as the union of all rate regions, C SDMA ({pfe}) defined in ©, over all values of pi and p 2 that satisfy 
XiPi + X2P2 = P*- Each C SD m({Pk}) is a pentagon with two vertices, each corresponding to one of the 
two possible decoding orders among the users [20]. The fact that user 2's message must be decoded 
first to achieve (Ri, R2) can be justified by the following observation. It is seen that the decoding order 
2 — ^ 1 always achieves higher rates for both users than the other decoding order 1 — ► 2, hence, it 
must be the optimal decoding order to achieve the rate target (R 1 , R 2 ) on the boundary of this capacity 
region. In general, for a deterministic SISO-MAC, all rate-tuples on the boundary of the capacity region 
under a weighted sum-power constraint can be achieved by a unique decoding order for each user. This 
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result is consistent with Theorem [T} i.e., given A^'s and h^s, the optimal decoding order of users can 
be first resolved by (1241) . and then their minimum powers can be found by (|23l ). 

Fading MIMO-MAC: Unfortunately, the contra-polymatroid structure is non-existent for the MAC 
when there is fading (e.g., the fading SISO-MAC) or there are multiple antennas at the receiver (e.g., 
the deterministic SIMO-MAC), and hence, it is non-existent for the general fading MIMO-MAC. 
Equivalently, the rate-tuples on the boundary of the capacity region under the weighted sum-power 
constraint for a fading MIMO-MAC, unlike the case of a deterministic SISO-MAC, might not correspond 
to a unique decoding order for each user. Fig. \5\ illustrates this fact by showing the capacity region of 
a two-user fading SIMO-MAC (ti — t% — 1, r — 2) under a sum-power constraint p\ + p 2 < 10 (i.e., 
Ai = A2 = 1). In this case, the channels {i2"i(z/)} and {H 2 {y)} are assumed to be independent vectors 
each having independent entries distributed as CA/"(0, 1). The capacity region for this case is shown to 
be symmetric. The dashed line and the dotted line show how two vertices of the constituting C SD ma({pa:}) 
sweep on the boundary of the capacity region as p 1 and p 2 vary while their sum is kept equal to 10. It 
is observed that the boundary rate-tuples of this capacity region indeed correspond to different decoding 
orders; e.g., point A corresponds to the decoding order 1^2, while point D does for 2 — >• 1. There is 
also a part of the capacity region (e.g., point B is in this region) that does not consist of any vertices. 
This part of the region is referred to as the time-sharing region and consists of the -45 degree boundary 
lines of the constituting Csdma({pa:})- Hence, any point in the time-sharing region is not achievable by 
successive decoding given any fixed decoding orders, and time-sharing the transmission rates and the 
decoding orders among two users is required. As a result, unlike the deterministic SISO-MAC, given 
Afc's and H k (u)'s, the optimal decoding order for the W-SPmin problem can not be resolved directly 
for this channel. 

One heuristic approach (e.g., [21]-[23]) for solving W-SPmin problem for the fading MIMO-MAC 
under SDMA might be searching through all possible K\ decoding orders and then finding the optimal 
decoding order that gives the minimum weighted sum-power to support the target rates. This approach 
might work for some special cases, but as explained below, is problematic in general. 
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For any fixed decoding order it, Problem \T\ can be written as 



A 



Minimize > Ai.Tr (Si.) 



(25) 



Subject to E^ — log 



1 Eti H ^(v) S k® H 1<M + 1 



(26) 



S k y VA; 



(27) 



where the left-hand-side (LHS) of (1261) is the achievable rate for MT ir(k), denoted as r^L, under the 
decoding order n. Except for MT 7r(l), the rate constraints in (1261) are not convex, rendering the above 
optimization problem non-convex in general and, hence, it can not be solved efficiently. A suboptimal 
method that approximately solves this problem is described as follows. Starting from the last decoded 
MT 7r(l), the method minimizes the power required to maintain the target rate for that MT, while 
considering MTs that have yet been decoded as interference. For example, for the two-user case and the 
decoding order of MT 2 followed by MT 1, the method first determines Si with the minimum p\ that 
satisfies Ej,[| log \H\{y)S\H\(y) + 1\] > R\ and then fixes Si and determines S2 with the minimum 
power p 2 that satisfies E^[| log \Hi (z^Siif^z/) + H 2 {p)S 2 H\{p) + I\] > R 1 + R 2 . Each of these 
two optimizations are convex and, hence, they both can be solved efficiently. The above algorithm is 
referred to as the greedy algorithm since each MT simply minimizes its own transmit power. For the 
special case of the fading SISO-MAC and SIMO-MAC where each Sk is simply a scalar and is equal 
to pk, the obtained p 1 and p 2 via the greedy algorithm are in fact optimal and minimize the weighted 
sum-power for any weights A under the given decoding order. This is because from (|26l) it can be 
shown that the minimum p n (k) required to support R n (k) for user 7r(/c) is an increasing function of the 
powers for the not-yet-decoded users, p w (i), . . . ,p n (k-i)- However, for the general fading MIMO-MAC, 
the above greedy algorithm might not be optimal because each MT now can adjust its covariance matrix 
to balance between minimizing its own transmit power and reducing the interference it causes to the 
users decoded earlier in the order. 

Nevertheless, even if the W-SPmin problem can be solved for each decoding order, the obtained 
powers that have the minimum weighted-sum among all decoding orders might still be suboptimal. 
This can occur when the target rate-tuple does not correspond to a unique optimal decoding order, e.g., 
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the rate-pair B in Fig. [5] that is on the boundary of the time-sharing region. 

From the above discussions, it follows that for the fading MIMO-MAC in general, the decoding order 
and the transmit covariance matrices for the MTs need to be jointly optimized for solving the W-SPmin 
problem under SDMA. This motivates the algorithm presented next. 

Proposed Algorithm: The proposed algorithm for Problem Q] is based on its Lagrangian [19], which 
is given below: 

K K 

£({S k }, {r k }, n) = Yl A * Tr (^) - - R ^ (28) 

k=l k=l 

where /x = . . . , fi K ) e denotes the vector of dual variables associated with the rate inequality 
constraints in ©. The variables, {S k } and {r k } belong to the set denoted by V, which is specified by 
the remaining constrains in © and (flOl) . Then the Lagrange dual function [19] is defined as 

= ,-, min ^({5 fc },{r fc },/x). (29) 

{S k ,r k }ev 

The dual function serves as a lower bound on the optimal value of the original (primal) problem, 
denoted by q*, i.e., q* > g(fx),Wn [19]. The dual problem is then defined as max^^ g(/x) [19]. Let the 
optimal value of the dual problem be denoted by d* that is achievable by the optimal dual variables /x*, 
i.e., d* = g(/J,*). For a convex optimization problem, the Slater's condition states that the duality gap, 
q* — d* > 0, is indeed zero if the primal problem has a feasible solution set [19]. By using sufficiently 
large user powers, the set Csdma ({£&}) in © can be made arbitrarily large to contain any finite rate 
target R as an interior point. In other words, we can always find a feasible set {S k } that satisfies any 
given rate constraint R for Problem [T] Thus, the Slater's condition holds and the duality gap is zero for 
Problem [TJ This result suggests that q* can be obtained by first minimizing the Lagrangian C to obtain 
the dual function g(fj.) for some given fi, and then maximizing g(fji) over all possible values ix. 

First, consider the minimization of C to obtain the dual function <?(//). In this case, /x is fixed and 
the variables are {S k } and {r k }. From (1281) . it follows that the minimization of C can be rewritten as 
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the following equivalent problem: 



Minimize 

{r k },{S k } 

Subject to 



K K 
k=l 



k=l 



k=l 

r E C smA ({S k }) (31) 
S k y VJfe. (32) 

From the definition of Csdma in ®, there are 2 K — 1 rate constraints implied by ((311) . which are difficult 
to be incorporated directly into the optimization. In order to simplify the problem, the following theorem 
is utilized to remove these constraints in (13TI) : 

Theorem 2: For any permutation n over {1, . . . , K} and a fixed set of covariance matrices {Sk}, 
defined as 



r% = K 



log 



Ei=l H n(i 


(z/)5 7r ( i 








)(f)5 ff ( 


J 7T(' 


)(») + * 



(33) 



is a verto: of the polymatroid C smA ({S k }) in M/f. Furthermore, for any p y 0, the solution to the 
following W-SRmax problem: 



K 

Maximize > p^rt. 

«=i 



Subject to r G C S dma({S7c}), 



(34) 
(35) 



is attained by a vertex r^*\ where 7r* is such that p n *(i) > P-k*{2) > • • • > P-k*{k)- 

Proof: Please refer to [13, Lemma 3.10]. ■ 
Notice that Theorem |2] holds for any given {Sk}- Furthermore, since minimization of — Y^ k p k r k is 
equivalent to maximization of 2~2 k p k r k , using Theorem [2j the rate constraints in (1311) can be safely 
removed and the problem in (1301) can be simplified as 



K K 

Minimize ^ A fe Tr(5 A; ) - ^ /M*0 E " 
fc = i fc = i 

Subject to S k y V/c, 



- log 

2 6 



(36) 
(37) 



where n is a permutation such that > • • • > p n (K)- By rearranging the terms regarding user rates 



in (1361) . the above problem becomes the minimization of 



K 



K 



/] AfcTr(5fc) — (/^r(fc) — ^7r(fe+l)) E v 



fc=l 



fc=l 



- log 

2 6 



\ H^){v)S n( i)Hl {i) (u) 

i=i 



(38) 
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with only optimization variables S k h 0, V/c, and n % nc+i) — 0. Since the above problem is convex 
with a twice differentiable objective function and positive semi-definite constraints, it can be solved 
numerically, e.g., by the interior-point method [19]. 

Next, the dual function g(/j,) is maximized over all possible values /x. The search of /x towards its 
optimal value /x* can be done, e.g., by the ellipsoid method [24], which utilizes the fact that the vector 
6, defined as 9 k = R k — r' h for k — 1, . . . , K, is a sub-gradient of g(fjt) at any n, where {S' k } and {r' k } 
are the minimizers of £({S k }, {r fc },/x) obtained via solving (|38l) . i.e., £({S' k }, {r' k },n) = g(fjt). 

Remark 4.1: It is noted that the algorithm proposed in [13, Algorithm 5.3] can also be modified to 
determine /x* for the problem at hand. However, from programming implementations, it is observed that 
this method may exhibit oscillation when some of /x£'s happen to be equal. In contrast, the ellipsoid 
method is more suitable because of its robust and superior convergence behavior. 

The complete algorithm for Problem [T] in the case of SDMA is summarized below. 

Algorithm 1: 

• Given an ellipsoid £[0] C M. K , centered at /x[0] and containing the optimal dual solution fx*. 
. Set i = 0. 

. Repeat 

1. For given fj,[i], solve the problem given in ([38]) to obtain an optimal solution set {S^fi]} and 
{ r fc[^]} that minimizes £({S k }, {r k },fi[i]) over V; 

2. Update the ellipsoid based on £[i] and the sub-gradients 9 k [i] = R k —r k [i], k = 1, . . . , K. 
Set fj,[i + 1] as the center of the new ellipsoid £[i + l]z 

3. Set % *- i + 1. 

• Until the stopping criteria for the ellipsoid method is met. 

One possible method to obtain the initial ellipsoid £[0] that contains the optimal dual solution fi* is 
given in Appendix U 

The primal-dual approach used for solving Problem [T] can be explained by the connection between 
the power region and the capacity region as described in Section [Till In Fig. [3] it is observed that 

2 Notice that when locating the center of a new ellipsoid, we need to add the constraint that (i y 0. 
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each power-tuple on the boundary of the power region for a given target rate-tuple, as the solution to 
Problem Q] for a given weight vector A, defines a capacity region that contains the target rate-tuple on its 
boundary. Consequently, the target rate-tuple can be expressed as the solution to a W-SRmax problem 
for some unknown user weight vector p over this capacity region. Clearly, the primal-dual approach 
establishes the above connection by finding the optimal dual variable /x* that is simply one candidate 
for the unknown weight vector p in this W-SRmax problem. From (1381) . it follows that the optimal 
decoding order of each MT is given by the magnitude of [i* k , i.e., the optimal decoding order n satisfies 
— ' " " ^ lKr(K) ■ Hence, the proposed algorithm successfully jointly optimizes the decoding order 
and the transmit covariance matrices of the MTs by exploiting the duality between the power region 
and the capacity region. 

Uniqueness of Solutions: So far, the proposed algorithm determines the optimal value of the primal 
problem q* (equal to that of the dual problem d*), the corresponding primal variables {S* k } and {r^} 
(it is yet claimed that these primal variables are the primal optimal solutions), and the dual optimal 
solutions fi* that satisfy, 

K K 

q* = d* = Y, A fc Tr(5*) - $>;(rj - R k ). (39) 

k=l k=l 

In the following, the issue on the uniqueness of these solutions is addressed. While uniqueness of the 
dual optimal variables fi* is not an issue in the convergence of the proposed algorithmic uniqueness of 
the primal variables, {S* k } and {r£}, plays an important role in obtaining valid solutions for the W- 
SPmin problem. Since a primal-dual approach is used, the obtained primal variables that minimize 
the Lagrangian at /x* might not satisfy the rate constraints in <£Sj) . Notice that these variables are 
minimizers of the Lagrangian and are not necessarily the primal optimal solutions. However, according 
to Karush-Kuhn-Tucker (KKT) optimality conditions [19], the primal optimal solutions also minimize 
the Lagrangian at /x*. Hence, if these Lagrangian minimizers can be proven to be unique, it follows 
that they satisfy the rate constraints in ([8]) automatically. 

3 By the primal-dual approach and the connection between the power region and the capacity region, it follows that /x* can be viewed as 
the weight vector p that attains the given rate requirements R as the solution to the W-SRmax problem over the corresponding capacity 
region. For the fading MIMO-MAC, as shown in Fig. [5] in general there is no "sharp" vertex with multiple tangent lines on the boundary 
of the capacity region under a weighted sum-power constraint. As a result, it can be inferred that /x* is unique for any rate-tuple on the 
boundary and, hence, the uniqueness of /x* is in general ensured. 
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Theorem 3: The primal optimal solutions {S* k } for Problem \T\ under SDMA is unique. 
Proof: Please refer to Appendix lU ■ 

If all /i^'s are positive and distinct, r* (e.g., shown by Point A in Fig. [5]) that maximizes J2kf x t r k 
over Csdma ({S^}) wml be one unique vertex of Csdma which itself is also unique according to 
Theorem [3] In this case, from the KKT conditions, r£'s automatically satisfy the rate constraints in ([8]). 
However, if /x£'s in some subset J C {1, . . . , K} are positive and equal, r^'s for the users in the set J 
may not be unique and consequently they may not satisfy the rate constraints in ([8]). This can be shown 
by point B, C and D in Fig. [5] where /i* = n%. Any point on the straight line between point C and D 
maximizes YliklAfk over C S dma because all these rate-pairs have the same sum-rate. Hence, if 
point B is the rate demand for our problem, because the simplified optimization in (1381) always tends 
to use a vertex of the capacity region as the solution for r*, the proposed algorithm would converge 
to either point C or D as r*, which clearly does not satisfy the rate constraints. However, this is not 
an issue in the convergence of the proposed algorithm. As far as {S* k } is unique, the target rate-pair is 
ensued to be some convex combination of at most K vertices of the unique C smk ({S* k }) . 

B. Solutions to Sum-Power Minimization Under Power-Profile Constraint 

This subsection presents the solutions to Problem [3l the SPmin-PPC problem under SDMA, where 
the transmit power of each MT is regulated according to a given power-profile vector a. The proposed 
algorithm is also based on a Lagrange primal-dual approach. The Lagrangian of the primal problem in 
(flTT) can be written as 

K K 

C(P, {S k }, {r fc }, fi, S) = P + 4(Tr(5 fe ) - a k P) - ^2 Vk(r k - R k ), (40) 

fc=i fe=i 

where /x = (fj,i, . . . , G M^" and 6 = (Si, ...,6k) £ ^+ are dual variables associated with the 
inequality constraints in (fT8l) and (|2TT) . respectively. Let T denote the set of primal variables specified 
by the remaining constraints in (fT9l) , (T20l) and (T22l) . the Lagrange dual function can be then expressed 
as 

(K \ K K 

1 - y^hat + y^5fcTr(S fc ) - V\ fc (r fc - R k ). (41) 
k=i ) k=i k=i 
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From (HTI) . it is necessary that 1 — J2 k=1 ^kO-k > for the dual function to be bounded from below. 
In the case of 1 — Y^k=i &k&k = ^' tne °pti ma l P mat minimizes the Lagrangian over T can take any 
positive value; while in the case of 1 — \\^ k=l ^k&k > 0, the optimal P must be equal to zero. Thus, 
in both case, the term associated with P in (HTI) is indeed zero and, hence, can be removed from this 
point onward. The optimal value of P can be then obtained as 

P* = max gin, 8). (42) 

/^0,<^0,E£Li<5 fc a fc <l 

Similar algorithm like Algorithm [TJ for Problem [TJ can be readily developed for solving this problem. It 
can be shown that for the problem at hand, g(fi, S) has sub-gradients, 6 and £ for p and 6, respectively, 
which are defined as 9 k = R k — r' k and ( k = Tr(S" fc ), k = 1,...,K, where {S' k } and {r' k } are the 
Lagrangian minimizers that satisfy £({S' k }, {r' k }, p, S) = g(p, 8). 

The convergence of the algorithm for this problem and the uniqueness of the solutions are similar 
as Algorithm [TJ Let {SI} and d* denote the corresponding primal and dual optimal solutions. It is 
worth mentioning here that the obtained primal solution for this problem P* might not be equal to 



the user sum-power £ 



K 

fern 



y Tr(S* k ) since some of the power constraints in (|2TI) . Tr(Sfc) < a k P, may 



not be active in generalo Actually, P* = J2k=i ^k^ T (^*k) shown as follows. Since P* > 0, it can be 
verified from (HTI) that J2k=i $k a k — 1- Moreover, by the KKT conditions for the power constraints, 
Ef =1 ^(Tr(S*) - a k P*) = 0. Hence, £?=i 5* k Tr(S* k ) = £f =1 5* k a k P* = P*. 

The primal-dual approach used for solving the SPmin-PPC problem is also based on the connection 
between the power region and the corresponding capacity region, similar as that for the W-SPmin 
problem. However, their difference lies in that for the SPmin-PPC problem, with reference to Fig. [31 
the weight vector for characterizing the solution on the boundary of the power region, A, and the 
weight vector for the corresponding capacity region, p, are both unknown and, hence, they need to 
be found under the given power-profile vector a as the corresponding optimal dual solution S* and 
H*, respectively. In contrast, for the W-SPmin problem, only the unknown p needs to be found as the 

4 For example, in the case of two-user deterministic SISO-MAC as shown in Fig. [4] (a), if the given power-profile vector a for Problem 
[3] is such that the intersected boundary power-tuple (aiP* , oli P*) of the power region is located on the vertical (or horizontal) boundary 
segment of the power region, the solutions to Problem [3] will converge to the upper (or lower) vertex (p*,P2) of the power region for 
which, clearly, pi = aiP* , but p* 2 < 0.2 P* ■ Thus, P* 7^ p\ + p 2 - 



20 

optimal dual solution fi* because A is already given. 

V. Characterization of Power Region for TDMA 

This section considers the characterization of the power region defined in © for the fading MIMO- 
MAC under TDMA, and for brevity only the W-SPmin problem (Problem [2]) is investigated. The 
alternative means for characterizing the power region based on the power-profile vector, i.e., the SPmin- 
PPC problem under TDMA, is omitted since it can be readily obtained given the techniques developed 
in Section IIV-BI for the case of SDMA. The proposed algorithm for Problem [2] is also based on the 
Lagrange primal-dual approach. 

First, it is noted that the constraints in (fT2l) and (fT3l) in Problem [2] can be combined and thus 
simplified, given the fact that in ((!]), the inequality constraints are always satisfied with equalities for 
power minimization. Hence, the Lagrangian of the primal problem can be written as 

K K 









jlog 





H k {u)^H\{u) + I 

Tk 



Ri 



(43) 



k=l k=l 

where /x = . . . , fx K ) G IR+ denotes the dual variables associated with the constraints in (fP2l) . The 
variables, {S k } and {r k }, belong to the set denoted by Q, which is specified by the remaining constrains 
in (fl4l) . (fT5l) and (fT6l) . Also note that the Lagrangian is a convex function of both {S k } and {r k }. By 
changing the variables as 



S k 



Tk 



For k 



K 

k=l 

K, define 



W k , k = 1, . . . , K, the Lagrangian can be rewritten as 

A' 



k=l 







(r k E u 


jlog 





H k (u)W k H\{u) + J 



— Ri 



(44) 



z k (W k ,fi k ) = \ k Tr(W k ) - ji k E v 



log 



H k (u)W k H[{u) 



(45) 



and 



z k (n k ) = min z k (W k ,Hk)- 
From (l44"l) . (|45T) . and (l46l) . the Lagrange dual function can be expressed as 



(46) 



{W k ,r k }eg 

K K 

= min ) T k z k (/i k ) + fi k R k 

{T k }-T k >0 Vk,J2k=i T fc =1 k= i k= i 



(47) 
(48) 
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The optimal value of the primal problem, denoted by q*, can be then obtained as 

Q* = max0(/j) (49) 

K K 

= E r **^)+E^- (50) 

k=l k=l 

Theorem 4: If {S* k }, {r k }, {r k } are the optimal primal solutions and are the optimal dual 

solutions for Problem [2] under the strictly positive weight vector A and rate target R, they must satisfy 

\kPt-nlr* k = <?Tl k=l,...,K, (51) 



where f k = Tr(S*), r* = r*E u \ log 



H k (v)^H[(v) + I 



Rk, and c* is a constant. 



Proof: Since R and A are both strictly positive, then so are the obtained solutions {r k }. As a result, 
it is necessary to have Zi(jt/J) = ■ ■ ■ = zk{[i*k), otherwise the minimization in (|48T) must lead to only 
one user assigned with the total time slot, i.e., T k , = 1 and r k = 0, k ^ k' where k' = argmin^ Zk(fJ>%). 
Using this equality and also (1451) . (l46l) . the proof is completed. ■ 

Using Theorem SI the algorithm for Problem [2] can be obtained as follows: In each iteration, the 
algorithm updates the dual variables such that Zi(jWi) = ■ ■ • = zr(hk) = c. It then checks whether 
the obtained rates can support more than the target rates, and increases c if they do or decreases it 
otherwise in the next iteration, until the rate targets are exactly met and c converges to c*. The details 
for the proposed algorithm are presented below. 

Algorithm 2: 



. Given c min < c* < c max . 
. Repeat 



!• C 2 (^min C ma x) 



2. For each k, obtain the optimal solutions W' k and fi' k such that z k (n' k ) = z k (W' k , n' k ) = cU Do 
the above for k = 1, . . . , K. 



3. Compute r' k such that R k = 7^E„ \ log H k (v)W' k H\{v) + / 
4 - If YLk=\ T 'k < !' c min c; otherwise c max <- c. 



for jfe = 



5 For each given /ik, Sfc(Mfc) can be obtained by minimizing ^(Wfc, /Ufc) over W k as in ((46} by means of a convex optimization 
method, e.g., the interior-point method [19]. The quantity fi' k for which Zk(fi' k ) = c can then be obtained by a bisection search over fit 
using the fact that Zk(fJ,k) is a decreasing function of fj,t- 
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• Until c max — c m i n < 5 where 5 is a small positive constant that controls the algorithm accuracy. 
Since Zk(fJ>%) < 0,Vfc, it follows that c* < 0. Thus, we can take c max = 0. Similar as Appendix HI we 
can obtain the upper bounds /4 s on the optimal dual solutions /i£'s. From (|43l) and using the fact 
that g(fi*) > 0, it is easy to show that c* > — Y^k=i l^t^-k- Thus, we can take c min = — ^2 k=i ^ Rk- 
At last, the convergence of the above algorithm as well as the uniqueness of the obtained solutions are 
ensured by the uniqueness of c* in Theorem HI 

VI. Numerical Results 

This section presents the power region for a fading MIMO-MAC with r = 2 receive antennas and 
K = 2 MTs each equipped with t\ = t 2 = 2 transmit antennas. It is assumed that the receive antennas 
at the BS are sufficiently separated that they experience independent fading, while the fading levels are 
correlated across the transmit antennas because of their realistic size limitations. Under this assumption, 

1 /2 

the employed channel model for MT k is given by H k {u) = H w {u)Q k for k = 1,2, where Q k G 
0kXtk denotes the transmit antenna correlation matrix for MT k and is assumed to be constant over 
all fading states of v. H w (y) E C rxtk denotes the Rayleigh-fading channel matrix that is independent 
across two MTs and across all fading states, and has independent entries distributed as £/V(0, 1). Similar 
as the proof given in [11] and [12], it can be shown that the expressions in (1381) and (|45l) . for SDMA 
and TDM A, respectively, are both minimized when the transmit signal covariance matrix, has the 
same set of eigenvectors as Q k , k = 1,2, i.e., if A k A k A\ is the eigenvalue decomposition of Q k , the 
optimal S k then takes the form of A k Y> k A\ for some diagonal matrix S fc . This observation reduces the 
number of (real) variables from t\ in S k to t k in T, k for MT k and in turn reduces the total algorithm 
complexityjj Monte-Carlo simulation with 5000 independent realizations of the random channels is used 
to approximate the actual expectation over the fading states. The simulation assumes that the target rate 
is R — [2 1] T nats/sec/Hz for two MTs, and 



Qi 



1 0.4 
0.4 1 



Q 2 



1 0.5 
0.5 1 



6 The proposed algorithms work for all kinds of CDI, e.g., with constant channel mean matrix, constant channel covariance matrix, or 
combinations of them in all general forms. However, for most of these cases, a similar variable-reduction like in this numerical example 
is not possible. 
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Fig. |6] and Fig. [7] show the power region for this fading MAC under SDMA and TDMA obtained by 
solving Problem \T\ and Problem [2l respectively. 

For the SDMA case as is shown in Fig. [6l there are two corner points on the boundary of the power 
region, denoted by point A and B, which can be obtained by the greedy algorithm described in Section 
HV-Al Recall that in this greedy algorithm, the BS first picks one possible decoding order for the MTs, 
and then starting from the last decoded MT, it minimizes the power required to maintain the target rate 
for each MT while considering the MTs that have yet been decoded as interference. In this figure, the 
power-pair A corresponds to the decoding order 2 — ► 1, while the power-pair B does for the reversed 
decoding order. The greedy algorithm achieves the optimal power-pairs for the W-SPmin problem under 
some weight vectors, e.g., point A for Ai 3> A2 and point B for A2 ^> Ai, but might be suboptimal 
for weight vectors other than these extreme choices. For example, the minimum value of OApi + 0.6p 2 
to support the target rate [2 1] T is 11.5 while the greedy algorithm leads to 12.8 and 13.3 units of 
power for power-pair A and B, respectively. Moreover, the boundary curve of the power region is not 
attainable by simply time-sharing these two corner points, as shown by the dashed line in Figj6l 

The power region under TDMA for this MAC is shown in Fig. [7J For comparison, the power region 
under SDMA is also included in this figure. This figure uses the log scale for the powers and, hence, 
the power region under SDMA looks as a non-convex set. The power savings achieved by SDMA are 
observed to be substantial compared to TDMA, even in this case where the number of transmit antennas 
at each MT is equal to that of the receive antennas at the BS, i.e., both SDMA and TDMA have the 
same number of degrees of transmission freedom in the spatial domain, which is two in this case. 
The power-pairs on the boundary of the power region under TDMA correspond to different time-slot 
durations, r k , k = 1,2, assigned to each MT. The power-pair A shown in Fig. [7]is achieved by assigning 
equal duration of time slot for both MTs, i.e., t\ = r 2 = 0.5 as in the conventional TDMA. Clearly, this 
power-pair is optimal for the W-SPmin problem under a unique weight vector A, and for the SPmin-PPC 
problem under a unique power-profile vector a, but it is suboptimal in all the other cases. For example, 
the minimum sum-power p\ + p 2 to achieve the target rate [2 1] T is 42 units of power for t\ = 0.66 
and r 2 = 0.34, as compared to 69 units of power for T\ = r 2 = 0.5. 
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VII. Conclusions 

This paper characterizes the power region for the fading MIMO-MAC. Motivated by a general 
relationship between the power region and the corresponding capacity region, the Lagrange primal-dual 
approach is employed to characterize all pareto optimal power-tuples on the boundary of the power 
region. These optimal power-tuples provide different power tradeoff among the MTs and also ensure 
the fairness of power consumption among them. The algorithms developed in this paper can be used 
in the wireless cellular network for the BS to control the transmit powers from the MTs in the uplink 
transmission. Two multiple access techniques, namely, SDMA and TDMA, are considered in this paper. 
It is observed that substantial power savings can be obtained by using SDMA compared to TDMA. 
This observation provides an important information-theoretic guidance for practical system designs, 
i.e., if the complexity for implementing the optimal SDMA can be tailed for, an enormous capacity 
gain is still possible over the conventional TDMA-based network. The multiuser transmit-covariance 
feedback scheme studied in this paper optimizes the transmit covariance matrices of all the MTs based 
on their long-term CDI. Hence, this scheme reduces significantly the feedback complexity compared 
to other feedback schemes based on the instantaneous channel realizations. As a result, this scheme 
is practically suitable for wireless channels that exhibit some consistent long-term channel statistics. 
The results obtained in this paper can provide insightful guidelines to many applications in wireless 
networks including resource allocation, partial channel feedback, and multiuser space-time code design. 

Appendix I 
Initial Ellipsoid for Algorithm Q] 

In the appendix, one possible method to obtain the initial ellipsoid S[0] for Algorithm Q] is presented. 
First, we obtain an upper bound ^ on fi* for any given j 6 {1, ... , K}. Let {S^} E V be any set 
of transmit covariance matricies that achieve {r^ } given by 

m ( R* k^j 
r? = < (52) 
{Rl + lk = j. 

From the definition of the dual function given by (|29l) , we have 

A 



g(S) < £({Sj°}, {rjp}, /O = A * Tr ffi) - /V 

k=l 
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Since g(fi*) > 0, it follows that 



K 



(53) 



k=l 



Thus, fj,f = ELi A * Tr ( S k ' )■ Similar upper bounds can be found for all other j. Next, S[0] can be 



chosen to cover the hyper-cube in M K specified by // 's. 

Appendix II 
Proof of Theorem [3] 

This Appendix proves the uniqueness of the solutions for the optimal transmit covariance matrices 
{S* k } in Problem [TJ Without loss of generality, it is assumed that /U* > - • • > fi* K > ij* k+ i = 0. If 
} and {S^} are two sets of optimal solutions for Problem [0 from (|28l) and by using Theorem 
[21 it follows that 



A' 



K 



K 



fc=i 



fc=i 



fc=i 



i=l 



(54) 



for j = 1,2. Since the problem at hand is convex, for any /3 G [0, 1], S k = fiS^ + ftS^ is also an 
optimal solution to satisfy (|54|) . where j3 = 1 — (3. This fact together with the concavity of the log | • | 
function implies that 



log 



(3A {1 \v) + {3A w (v 



(2), 



-/?log 



-/Slog 



0. 



(55) 



where A U) (u) = £f =1 H k (u)S k j) H k {u)^ + i". Let /(/3) denote the function on the LHS of the above 
equation, then f((3) — 0, for all < @ < 1. Because f{(3) is twice continuously differentiable, both of 
its first and second derivatives must vanish, i.e., 

d 2 f{(3) 



d(3 2 



Tr (v) - A {2) (u) pA w (v) + @A 



v 



0, Wf3. 



(56) 



For every v, the matrix in Tr(-) of the above equation is a positive semi-definite matrix and, hence, it 
has a nonnegative trace. Since the expectation of a nonnegative random variable is zero, it must be zero 



i.s., or A^\u) = A^ 2 \i>) a.s., which implies that S 



(i) 

k 
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Fig. 1. The fading MIMO-MAC with the multiuser transmit-covariance feedback. 




Pi 

Fig. 2. Characterization of power region via WSPmin or power-profile vector. 
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Fig. 3. The relationship between each boundary point of the power region and its corresponding capacity region for a two-user fading 
MAC. (a) Power region under the rate constraint (Ri,R,2), where (pi,pl), represented by point A, is on its boundary and achieves the 
minimum weighted sum-power under the weights Ai and A2, i.e., Xipl + A2P2 = P""\ (b) Corresponding capacity region of (p\,pi) under 
the weighted sum-power constraint Aipi + A2P2 < p*, where (Ri, R2), represented by point B, is on its boundary. Moreover, (Ri, R2) 
maximizes the weighted sum-rate pin + p2»"2 for some nonnegative weight vector p. 




Fig. 4. The relationship between each boundary point of the power region and its corresponding capacity region for a two-user 
deterministic SISO-MAC under SDMA. (a) Power region under the rate constraint (.Ri,i?2), where (p*,p2) is on its boundary and 
achieves the minimum weighted sum-power for the weights Ai and A2, and the associated decoding order is 2 — > 1 since j± > j^; (b) 
Corresponding capacity region of (pi, pi) under the weighted sum-power constraint Aipi + A2P2 < p*, where p* = Xipl + \2pl- 
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Fig. 5. Capacity region for a two-user symmetric-fading SIMO-MAC under a sum-power constraint: pi + p-2 < 10. The multiple access 
technique is SDMA, and ti = ti = 1, r = 2. 
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6. Power region for a two-user transmit-correlated fading MIMO-MAC under SDMA with ti = ti = 2, r = 2, and the target rate, 
[2 1] T nats/sec/Hz. 
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Fig. 7. Power region for a two-user transmit-correlated fading MIMO-MAC under SDMA and TDMA with ti = ti = 2, r = 2, 
the target rate, R = [2 1] T nats/sec/Hz. 



